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Kinematics of hypervelocity stars in the triaxial halo of the 
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ABSTRACT 

Hypervelocity stars (HVSs) ejected by the massive black hole at the Galactic center 
have unique kinematic properties compared to other halo stars. Their trajectories 
will deviate from being exactly radial because of the asymmetry of the Milky Way 
potential produced by the flattened disk and the triaxial dark matter halo, causing 
a change of angular momentum that can be much larger than the initial small value 
at injection. We study the kinematics of HVSs and propose an estimator of dark halo 
triaxiality that is determined only by instantaneous position and velocity vectors of 
HVSs at large Galactocentric distances (r > 50 kpc). We show that, in the case of a 
substantially triaxial halo, the distribution of deflection angles (the angle between the 
stellar position and velocity vector) for HVSs on bound orbits is spread uniformly over 
the range 10°— 180°. Future astrometric and deep wide-field surveys should measure the 
positions and velocities of a significant number of HVSs, and provide useful constraints 
on the shape of the Galactic dark matter halo. 
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1 INTRODUCTION 

Recent observations have revealed the existence of a popu- 
lation of hypervelocity stars (HVSs) traveling in the halo of 
the Milky Way (MW) with Galactic rest-frame velocities v r f 
in the range between +400 and +750 km s 



Brown et al 



Brown et al 



| 2005l:|Edelmann et al 1 l2005l : lHirsch" et al. 2005; 
l2006al lbl). HVSs are probably B-type main sequence stars 
with lifetimes < 100 Myr, Galactocentric distances > 50 
kpc, and move with speeds large enough to escape from the 
Galaxy. The significant excess of B-type stars with velocities 
+ 275 < f r f < +450 km s -1 and distances > 10 kpc observed 
bv lBrown etafl (|2007h may also be an indication that many 

HVSs are ejected into the halo on bou nd orb its. 

HVSs were first recognized by iHillsl (|l988l ) as an 
unavoidable byproduct of the presence a massive black 
hole (BH) in the Galactic center. Only a close encounter 
with a relativistic potential well can accelerate a 3-4 Mq 
star to such extreme velocities, and at least three dif- 
ferent ejection mechanisms have been proposed: the in- 
teraction between background stars and an intermediate- 
mass black hole (IMBH) inspiralling towards Sgr A* 



* Also a Hubble Fellow at the Department of Astronomy, Uni- 
versity of California at Berkeley, Berkeley, CA 94720. 
f Also at the Max-Planck-Institut fur Astrophysik, Karl- 
Schwarzschild-Str. 1, 85740 Garching, Germany. 



jYu fc Tremainell2003l ; lLevinll2006l : iBaumgardt et alj|2006l : 
ISesana et all I2006T ). the disruption of stellar binaries in 
the tidal field of Sgr A* dHillsl Il988l: lYu fc Tremaind 120031 : 
iGinsburg fc Loebll2006l ; iBromlev et all 120061 ) . and the scat- 
tering o f stars off a cluster o f stellar-mass BHs orbiting 
Sgr A* jO'Learv fc Loebll2006l ) In all these models, HVSs 
have unique kinematics compared to other halo stars: 1) 
they have almost zero initial specific angular momentum at 
ejection, ^/GM m r p ~ 4.0 x 10" 6 kpc 2 Myr _1 (MB H /3.6 x 
lO 6 M ) 1/2 (r p /lO~ 6 kpc) 1/2 , where Mbh the mass of Sgr 
A* and r p the pericenter distance of the star; 2) their high 
speeds diminish the impact of two-body relaxation or dy- 
namical friction effects on their motion; and 3) their tra- 
jectories will deviate from being exactly radial because of 
the asymmetry of the Milky Way potential produced by the 
flattened disk and the triaxial dark matter (DM) halo, caus- 
ing a change of angular momentum that can be much larger 
than the initial small value. (For reference, a 1 kms^ 1 devi- 
ation of the velocity from the radial direction at 50 kpc rep- 
resents a change of 0.05 kpc 2 Myr -1 in specific angular mo- 
mentum.) Proper-motion measurements of HVSs may there- 
fore become a key diagnos tic tool for constrai ning the shape 
of the Galactic potential dGnedin et al.ll2005l ). 

Triaxial halos are a generic prediction of the hierarchi- 
cal, cold dark matter (CDM) models of structure forma- 
tion. Dissipationless cosmological simulations typically pre- 
dict minor-to-major density axis ratios in the range 0.4-0.8 
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(e.g. IJing fc Sutol 12002 V with the asphericity of the po- 
tential increasing ra pidly towards the center of the halo 
l|Havashi et al. | [2006h . Gas cooling tends to circular ize the 



potential (e.g. Dubinskill994l ; lKazantzidis et al. 1120041 ), while 
subsequent merge rs produce highly elongated remnants (e.g. 
M oore et ai1l2004h . Studies of weak gravitational lensing and 
X-ray observations of elliptical galaxies show that halos are 
significantly flattened, i n fair agreement with results from 
nume rical simulations l|Hoekstra et al.l |2004| ; iBuote et al.l 
2002). Yet the coherence of tidal debris from the Sagittarius 
dwarf galaxy appears to indicate that the inner halo of the 
MW is nearly spherical and therefore in conflict with CDM 
predictions (Ibata et al. 2001; but see Helmi 2004). 

In this paper, we study the kinematics of HVSs in the 
MW as a probe of the triaxiality of the Galactic halo. The 
outline is as follows. In §[2l we analyze the motion of HVSs in 
a flattened or triaxial gravitational potential. We provide a 
concise statistical estimator for the triaxiality of the Galactic 
halo potential through the measured angular momenta of 
HVSs. In § [21 we review the Galactic potential model to 
be used in our calculations. In § U we perform numerical 
simulations of the motion of HVSs to study their kinematics. 
Finally, in §[S] we summarize our conclusions. 



2 MOTION OF HYPERVELOCITY STARS 

Consider a star with position vector r moving with velocity 
v in a gravitational potential $ = $ S ph(>") -1- & nS ph(x, y, z), 
where $ sp h and $ nsp h are the spherically-symmetric and 
aspherical component of the the potential, (x,y,z) are 
Cartesian coordinates, and r — \/ x 2 + y 2 + z 2 . The rate 
of change of the specific angular momentum of the star, 
J = r x v, is equal to the torque, 



dJ/dt = - fx V<3> = — f X V<& nS ph, 
and has components 

dJ x /dt = 2yz (d$ nBp h/dy 2 - d$ nsph /dz 2 ) 



dJy/dt = 2xz (d&ns P h/dz 2 — d$ n s P h/dx 
dJ z /dt = 2xy (d$ na p h /dx 2 - d$ nsp h/dy 



(1) 



(2) 



It is convenient to change from Cartesian to spherical co- 
ordinates, (x,y,z) — (r sin 9 cos <f), r sin 6 sin cj>, r cos 9), and 
combine the above equations to yield 



dj x jdt 



sin U sm cos ( 



+ 



dJy/dt 



sm a cos cos ( 



+ 



dJz/dt 



sin a cos sin < 



= 0. (3) 



From the definition of angular momentum it is also easy to 
derive 

Jx+Jy + Jz = 0, (4) 

where 

{Jx = J x j (sin 9 sin <f> cos 9) , 
J y = J y / (sin 6 cos (f> cos 9) , (5) 
Jz = J z / (sin 2 9 cos 4> sin (f>) 

are determined directly from the position and velocity of 
the star. Note that equations ((3]) and Q are rotationally 
invariant, that is, they do not change when arbitrary rota- 
tions are applied to their arguments. Below we apply the 
above analysis to the motion of stars in two simple cases of 
non-spherical potentials. 



• If the non-spherical component of the gravitational po- 
tential is axisymmetric about the plane 2 = 0, 



$nsph = $nsph(.R = \A 2 +y 2 , Z), (6) 

then <9<I> nS ph(-R, z)/dx 2 — <9<l? nsp h(-R, z)/dy 2 , and J z is con- 
served. Stars ejected from the Galactic center on radial or- 
bits move in a plane with 



Jz = 0. 



(7) 



• If the non-spherical component of the potential is tri- 
axial, 

$ns P h = $ nsp h(z 2 + y 2 /p 2 + z 2 /q 2 ), (8) 
then a triaxiality parameter can be defined as 

^^4- (9) 



If p = q = 1, the potential reduces to the spherical case. If 
p = 1 and q ^= 1 (T = 0), q = 1 and p^l, ot p = q ^ 1 
(T — 1), the potential is axisymmetric. If q < p < 1, the 
triaxiality parameter is < T < 1. In a triaxial potential, 
equation ([2j can be written as 



dJz/dt 



dJy/dt 



sin U cos sm < 



sin a cos cos f 



(10) 



For HVSs moving away from the Galactic center on radial 
orbits, the deviation of their trajectory from the initial ejec- 
tion direction, (59,Stp), is small. Replacing the angles (9,(f>) 
in equation (|10|l with (9 ± 59, <f> ± 8(j>) and integrating yields 



T = 



Jz_, 259 
~J y [T sin(20) 



tan < 



+ 



+ 



5 2 4> 



25950' 



sin(26') tan < 



+ .-], 



(11) 



where the sin and cos arguments have been kept con- 
stant in the integration. The term in parenthesis specifies 
the maximum systematic error on the triaxiality estimator 
T = —Jz/Jy, and numerical calculations (see § 4 below) of 
the motion of HVSs in a triaxial potential show that the typ- 
ical error on T is actually smaller. For a sample of iV HVSs 
we can use equation (TTTJ to construct a statistical estimator 
of triaxiality as 



T 



with standard deviation 

-1/2 

/v^ 1 \ 

CTrp = 




Here 



T, 



Jz,i 

Jy,i 



J y< i sin ( 



(12) 



(13) 



(14) 



and OTi is its error. Note that the parameters Ti are fully 
determined by the instantaneous positions and velocity vec- 
tors of the HVSs in the sample, and that, while the ra- 
tio — Jz,i/ Jx.i can also be used to estimate the triaxiality 
(1— p~ 2 )/(g~ 2 — p~ 2 ), this does not provide any independent 
information since it can be trivially derived from J z ,i/Jy,i 
using equation Q. 
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In the following we set the z-axis of the triaxial potential 
to be normal to the Galactic disk plane, and denote with 
rjo the angle measured counter-clockwise from a reference 
direction (e.g. the line from the Galactic center to the Sun) 
to the a;-axis of the potential. The ratio —JlilJLi in a 
frame forming an angle r\ with the reference direction can 
be written as 



7" 

3y,i 



where 



A = 



14- 



-1 B 



uin(2#) 



sin [2 (77 — 770 + > 



B = 



sin(j7 — 770) sin 4>'i 
cos(r? - r)o + 4>'i) 



(15) 



(16) 



If 77 7^ rjo, the values — J' z ,il J' y ,i ma y spread out over a large 
range due to the different angles in the sample. The angle 
r]o can be estimated by minimizing the weighted variance of 
T' = - 7' ■ 1.7! ■ ■ 



N 

1=1 



T- - T> 



(17) 



where a T > is the error of T[ and the values of T" and crj, are 
obtained from equations f| 1 2 p and (|13[) . 



3 GALACTIC GRAVITATIONAL POTENTIAL 

We use here a four-component model for the gravitational 
potential of the Milky Way, g> = ^BH + ^buigc + ^disk + ^haio, 
where fcf. iGnedin et aJjl200a ): 



$bh is the contribution of Sgr A* , 
GMbh 



(18) 



with mass M B h ^ 3.6 x 10 6 M |Ghez et alj 120051 ; 
lEisenhauer et al.ll2005l ). The radius of influence of Sgr A* is 
GM m /al ~ 1.6 pc (Mbh/3.6 x 10 6 M Q )(100kms- 1 /cr c ) 2 , 
where a c is the one-dimensional stellar velocity dispersion 
in the Galactic center. 

• ^buige is t he contribution of the spherical bulge 
l|Hernauisdll990D . 



GM] 



bulge 



bulge 



Q-bulgc 



(19) 



with mass Afbui ge = 10 10 M Q and core radius ab u i gc = 
0.6 kpc. As both the bulge mass and size are small compared 
to those of the Galactic disk and halo, a slight deviation of 
the bulge from sphericity will not have a significant effect 
on the change of angular momentum of HVSs in the halo. 

• *£disk is the contrib ution of the axisymmetric disk 
jMivamoto fe Nagailfl975l ). 



$ diak (R,z) 



GM dis 



(20) 



R 2 + ( ddisk + yjz 2 + b 2 diB 



with mass Afdisk = 4x10 Mq , scale length adisk = 5 kpc, 
and scale height 6disk = 0.3 kpc. 

• ^haio is the contribution of the triaxial dark matter 



halo, 

$hai (x,y,z) = $NF\v(r'), 



2 2 \ V2 

1/3 1/3/2.2/ Z 



<E > NPw(r' t ) 



GM 200 ln(l + r f /r s 
r a f(c) r*/r a 



(21) 
(22) 

(23) 



Here c = r2oo/f s is the halo concentration parameter, 7-200 
the radius within which the enclosed average density is 200 
times the mean matter density, r s the scale radius, and 
f(c) = ln(l + c) -c/ (l + c). This generalization of an NFW 
l|Navarro et all 1 19961 ) model ensures that the spherically- 
averaged potential of the triaxial halo is similar to that 
of a sp herical halo with t he same mass M200 and scale 
radius (jHavashi et al.l 120061 ) . We choose the following pa- 
rameters for the Milky Way halo: c = 15, r2oo = 389 kpc, 
M200 = 1.8 x 10 12 M (e.g. iDiemand et alJl2007t >. and as- 
sume that the z-axis of the halo and disk potentials coincide. 
We also assume, for simplicity, that the asphericity of the 
potential is constant with radius r, and set p = 0.8 and 
q = 0.7 in the calculations below. 

CDM halos are not smooth but have a wealth of 
substructure on all resol ved mass scales (e.g. iMoore et al.l 
Il999l : iKlvpin et al1ll999h . The assumption made above of 
a smooth gravitational potential neglects the deflection of 
HVS trajectories by halo substructure. A star with velocity 
v passing within a distance r* from a subhalo of mass M au h 
will change its velocity by an amount Sv = 2GM su ^/(vr t ). 
The mass function of substructure in a Milky Way-sized halo 
can be described as N(> M Buh ) = 6.4 x 10" 3 (Af 2 oo/M sub ) 
in the subhalo mass range 10 6 M < M sub < ,fma X M2oo 
ijDiemand et all 12007? ) . with / max = 0.01. The probability 
that a HVS ejected from the Galactic center is subject to a 
velocity deflection > Sv can then be estimated as 



P(> 5v) 



' 200 

7. 



dN 
dM su] 



V 0.01 / v 



lO^kms" 



-dM sut 



V c 



140 kms" 1 



v 



1 km. 1 



Sv 



(24) 



where V c = (GA^oo/^oo) 1 ^ 2 is the halo circular velocity. 
This probability is quite low so that a smooth potential is a 
good assumption in this work. 



4 NUMERICAL CALCULATIONS 

In this section we perform numerical calculations of the mo- 
tion of HVSs in the MW gravitational potential. 



4.1 Initial conditions 



According to the study of lYu fc Tremainel (|2003l ) . three- 
body interactions between ambient stars and a BH pair 
(where the second ary BH may be an IMB H inspiralling to- 
wards Sgr A*; e.g.. lBaumgardt et al1l2006T ) expel HVSs (v > 
10 3 kms _1 ) at a rate that can be as large as ~ 10 _4 yr _1 
(for a binary with semimajor axis 0.5 x 10 -3 pc and mass 
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Figure 1. Present-day spatial distribution and velocity vectors of HVSs ejected by a binary BH at the Galactic center. The length of 
each vector is proportional to speed. Ten thousand HVSs were generat ed at a constant rat e in the past 10 9 yr, with a velocity and spatial 
distribution obtained from the three-body scattering experiments of ISesana et al. | ll2006h . The binary has a mass 3.6 X 1O 6 M , mass 
ratio of 1/81, orbital semimajor axis equal to O.la^, and eccentricity 0.3, and orbits in the (x,y) plane. The reference axes are set along 
the (x, y, z)-axes of the triaxial halo potential in eq. (122[ i. 



ratio of 0.01). Tidal break-up of binary stars ("Hills' mech- 
anism") ejects HVSs at a rate ~ 10 — 5 (77/O.I) yr -1 , where 
r\ is the fraction of stars in binaries with semimajor axis 
< 0.3 AU. Close encounters of two single stars in the vicin- 
ity of Sgr A* may also produce HVSs but at the negligible 
rate of 10 _11 yr _1 . Different ejection mechanisms give ori- 
gin to different spatial and velocity distributions of HVSs. 
While Hills' mechanism predict HVSs to be expelled isotrop- 
ically at an approximately constant rate, in models involv- 
ing a BH pair HVSs are ejected preferentially within the 
orbit al plane of the binary in a short burst lasting a few 
Myr (| Zier fc BiermanrfeOOll ; lLevinll2006l ; ISesana et atlbOOrj. 
2007). In the latter case the degree of anisotropy depends 
on binary separation, the mass ratio, and the orbital eccen- 
tricity of the BH binary. 

For illustrative purposes, we assume in this paper that 
HVSs are ejected from the Galactic center by the BH bi- 



nary mechanism. We use the stellar spatial and velocity 
distributions deriv ed from the scattering experiments of 
ISesana et al.l {2006) for a binary with mass ratio M2/M1 = 
1/81, semimajor axis a — 0.1 ah, and eccentricity e = 0.3. 
The binary orbit is in the Galactic disk plane, and the 
velocity of the lighter hole at pericenter is direct ed along 
6 = 3n/2. The "hardening" radius ah is defined as l|Quinlanl 
Il99fj ) 



GM 2 




When a > ah, the binary separation decreases both by 
dynamical friction and three-body interactions with low- 
angular momentum stars passing in its immediate vicin- 
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Figure 2. Panels (a)-(c): Deflection angle a = arcsin(|rx v\/rv) between velocity and position vectors for all HVSs plotted in Fig.[T] as 
a function of Galactocentric distance r. The panels show different velocity ranges or different scales in a. Panel (d): transverse velocity in 
the Galactocentric frame versus distance r for all stars shown in Panel (c). Panels (a)-(d): triaxial halo potential withp = 0.8 and q = 0.7 
(see eg. 1221 1. Panels (e)-(h): spherical halo potential with p = q = 1. Different colors depict different velocity ranges: 300-400 kms 
(magenta), 400-600 kms" 1 (red), 600-SOOkms- 1 (blue), 800-L 000 km s" 1 (green), 1 000-1200 km s" 1 (cyan), > 1200kms- 1 (black). 
HVSs with angles between 5° and 180° in Panels (a), (b), (e) and (f) are bound stars with significantly bent orbits, and their detailed a 
distribution depends on halo triaxiality. Stars with small deflections (a < 5°) are either unbound (and have distances > 10 kpc) or are 
in the initial phases of their orbital period (and are closer to the Galactic center, see also Fig. [5}. The large cross in Panel (a) flags the 
locus of the bound star whose orbit is shown in Fig. [3] 
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Figure 3. The bound orbit of a simulated HVS ejected from 
the Galactic center. The orbit is significantly bent, rather than 
radial as expected in a spherical potential. The star's travel time 
820 Myr, and its present-day velocity and Galactoccntric distance 
are 300 kms" 1 and 10 kpc, respectively. 
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Figure 4. Histogram of the deflection angle distribution of all 
stars shown in Fig. [2j a ) (solid line) and (e) (dotted line) with 
r < 20 kpc and 300kms _1 < v < 400 km s" 1 . The bound stars 
are clustered around 0° and 180° in a spherical halo potential 
(dotted line), while their distribution is spread uniformly over the 
range 10°-180° in a triaxial halo potential (solid line). 



ity. After the binary becomes "hard" (a < a^), the 
bound pair loses orbital energy mainly through three- 
body interactions until gravitational radiati on takes over 
jBegelman. Blandford fc Reeslll980l ; IYu1|2002| ). The ejection 
speed of the stars at infinity spans a range with r.m.s. ~ 
7 x 10 2 kms" 1 (A/2/0.01 MBH) 1/2 (10~ 3 pc/a) 1/2 (M B H/3.6 x 
10 6 M ) 1/2 (see eq. 33 in lYu fc Tremaineil2003l ). In the cal- 
culations below, we assume 10 4 HVSs are ejected from the 
Galactic center at a constant rate over the last 10 9 years, 
and ignore for simplicity the orbital evolution of the binary 
during such timescale. Different stars move independently in 
the Galactic potential. A fourth-order Runge-Kutta method 
with adaptive stepsize control is used to solve the differential 
equations of the motion of the stars. Note that the adopted 
mass ratio and semimajor axis are within the allowed pa- 
ram eter space for a BH pa ir at the Galactic center (see Fig. 
2 in lYu fc Tremainell2003h . 



4.2 Results 

We use the initial conditions described above and numeri- 
cally integrate the orbits of HVSs in the Galactic potential. 



Figure 5. Deflection angle versus travel time from the Galactic 
center. Top panel: all stars shown in Fig. [2{a) with r < 50 kpc 
and 300 kms -1 < v < 400 kms -1 . Bottom panel: all stars shown 
in Fig. [2]» with r < 200 kpc and 400 kms -1 < v < 600 kms -1 . 
HVSs with short travel times, < 50 Myr (top) and < 400 Myr 
(bottom) have a angles smaller than 5°. The large cross in the 
top panel flags the locus of the bound star whose orbit is shown 
in Fig. [3] 



Figures [T] and [2] show maps of stellar position and veloc- 
ity vectors at the present time and their deflection angles 
a = arcsin(|r x v\/rv). For HVSs with v > 600 kms -1 , de- 
viations are quite small, a < 5°, at all distances within 
200 kpc. Lower velocity stars at small distances can in- 
stead be bound, and their deflection angles extend to 180°. 
Many of the stars with 300 kms -1 < v < 400 kms" 1 and 
r < 20 kpc in Figure [2ja) follow bound trajectories, while 
no stars in this velocity range have substantial deviation 
angles at large Galactocentric distances. One example of a 
bound orbit is shown in Figure [3] where the trajectory has 
been significantly bent by the triaxial halo and the flattened 
disk potentials, and the star does not return to the Galactic 
center. Note that bound stars with v > 400 kms -1 and large 
a's are typically observed at smaller distances (r < 10 kpc) 
than bound stars of lower velocities. This is because a larger 
ejection speed from the Galactic center is needed to main- 
tain a high velocity at large distances, and such stars may 
then have either escaped from the halo or not have had suf- 
ficient time to come back and show a significant bend in 
their orbits. (For reference, the local escape speed is about 
500 - 600 kms -1 . ISmith et all [200^ .1 We find that 60% of 
the stars with v > 300 kms" 1 at r < 20 kpc have velocities 
v < 400 kms" 1 , and 20% of the stars with v > 300 kms" 1 
at r < 10 kpc have velocities v < 600 kms -1 . 

Figure [5] shows that bound stars with large deflection 
angles (5° < a < 180°) have generally traveled a long 
time after ejection (>50 Myr for r < 20 kpc), and many 
of them have experienced at least one orbital period (Fig. 
[3|. HVSs with small deflections (a < 5°) have instead a 
short travel time and are generally on the initial phases 
of their first orbital periods (see the concentration of the 
stars at the left bottom of the panels in Figure [5}. The 
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Figure 6. The variables —JzlJy and (vfr)J'y versus Galactocentric distances r for all HVSs far from the plane of the disk, \9— 90° | < 30° . 
Stars near the plane are included only in Panel (d) . Different colors represent different velocity ranges as in Fig. f2] Panel (c) shows J y 
on a different scale than Panel (b), and the scattered magenta dots are bound stars with significantly bent orbits. Panel (e) shows the 
effect of removing the disk potential from the calculation. Panel (f) shows the results in a reference frame with r\ — rjo = 10°, different 
from the triaxial frame of the halo potential. The solid line represents the triaxiality of the halo potential T = 0.54 (with p = 0.8 and 
q = 0.7 in eq. [9} assumed in the calculation. 
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transverse velocities (in the Galactocentric frame) of stars 
with v > 300 km s -1 are typically higher than 3 kms -1 , and 
can be up to 30 kms" 1 (hundreds of kms" 1 ) for unbound 
(bound) HVSs. Note that 3 kms" 1 corresponds at a dis- 
tance of 100 kpc to a proper motion of 20 /ias in three years, 
which can be resolved by the next generation of astromet- 
ric surveys like GAIA. According to our calculations most 
HVSs would have transverse velocities > 3 kms -1 even in 
the case of a weakly triaxial halo with (p, q) — (0.95,0.9): 
these transverse velocities are larger than those associated 
with the bending of stellar trajectories caused by the ax- 
isymmetric disk (see Fig. 0r) . 

A comparison between Figures HJe)-(g) and |2ja)-(c) 
shows that the distribution of deflection angles versus dis- 
tance is different for trajectories in a spherical rather than 
triaxial halo. In Figure [2je) and (f), the a angles of bound 
stars are clustered around 0° and 180° because their orbits 
are highly eccentric with little bending (see also Fig. Q. 
By contrast, in Figure [2j a) and (b) the deflection angles of 
bound stars lies at intermediate values. This difference can 
be used as an indicator of the triaxiality of the MW dark 
matter halo. 

Figure EJa)-(c) depicts the values of — JzjJy and 
(v/r)J y versus distance of all the HVSs plotted in Figure 
[T] having v > 300 kms -1 and far from the plane of the disk, 
i.e. with |6>-90°| > 30°. Figure^d) shows all stars with 9 in 
the range 0°-180°. As seen in Figure[^a), the ratio —JzjJy 
at large distances (r > 100 kpc) is close to the triaxiality of 
the halo potential (T — 0.54, solid line) assumed in our cal- 
culations, with only a small scatter. The scatter is larger in 
the quantity (v/r)J y plotted in Figure H^b). The slight off- 
set (< 10%) of the stellar dots from the solid line is partly 
due to the disk potential that causes an additional bending 
of stellar trajectories towards the disk plane, increasing \J y \ 
without changing J z (see the Appendix for details about 
the correction of such offset owing to the Galactic disk). As 
shown in Figure |6][e), HVSs become better tracers of triaxi- 
ality after removing the disk potential from our calculations. 
The higher the velocities, the smaller the offset. The curva- 
ture of the ratio —JzjJy towards smaller values at small 
distances is also an effect of the disk. The ratio — Jz/Jy for 
stars close to disk plane is not a good approximation of halo 
triaxiality even at large distances (see Fig. [6jl) . We plot the 
values of J z versus J y in Figure [7] where Panel (a) shows 
stars with \9 — 90° | < 30°, and Panel (b) stars with 8 in 
the full range 0°-180°. The effect of the disk potential can 
be seen in the spreading of the stellar dots from the solid 
line Jz — —TJ y to the dashed line J z = 0. Note that the 
scatter in the ratio —Jz/Jy increases significantly if the ref- 
erence frame is different from the frame of the triaxial halo 
potential (Fig. [(J). 



5 SUMMARY 

We have studied the unique kinematics of HVSs ejected from 
the Galactic center with almost zero initial specific angular 
momentum. HVSs can travel in the Galactic halo on either 
bound or unbound orbits, and their spatial and velocity dis- 
tribution at large Galactocentric distances (r > 50 kpc) con- 
tain information on the asphericity of the halo gravitational 
potential. We have proposed an estimator of the triaxiality 



of the Galactic dark matter halo that is determined solely 
by instantaneous position and velocity vectors of HVSs, is 
independent of the details of the ejection mechanism, and 
does not require an accurate knowledge of halo mass. Future 
astrometric and deep wide-field surveys of HVSs should de- 
tect significant numbers of HVSs, which could be used to 
determine the triaxiality of the MW halo by applying the 
method proposed in this paper. 

The new class of possibly bound HVSs with veloc- 
ities +275 < v T { < +450 kms" 1 recently observed by 
iBrown et all (|2007T l has Galactocentric distances in the 
range 30-60 kpc or 10-20 kpc depending on whether they are 
main-sequence or blue horizontal branch stars. In the first 
case (main-sequence stars at large distances), they have a 
lifetime of < 100 Myr and are, according to Figure [SJ on the 
initial phases of their first orbital periods. Their deflection 
angles are expected to be rather small, supporting the fact 
that a significant excess of B-type stars is observed only at 
large positive velocities (|Brown et alj lioo?^. If bound HVSs 
are blue horizontal branch stars instead at smaller distances, 
the travel time of the stars can be much longer than 50 Myr 
as at the ejection moment the stars may not necessarily be 
blue horizontal branch stars but at some pre-blue-horizontal- 
branch stage. Such stars may have experienced at least one 
orbital period, and many of them should be returning to 
the Galactic center or their orbits should have been signif- 
icantly bent by the asymmetric Galactic potential (Fig. [SJ). 
This scenario appear unlikely since it does not agree with 
the observed positive radial velocities. 

It is interesting at this stage to provide an example of a 
statistical estimate of halo triaxiality from a mock sample. 
Let us assign a measurement error of a a = 3 kms -1 (a — 
x, y, z) to the one-dimensional velocities of all stars in Figure 
IT a) having v > 300 kms" 1 , |6> — 90° | > 30°, and 55 kpc < 
r < 200 kpc (our error analysis assumes that the distances to 
HVSs are known to within ten percent). We have simulated 
the observed velocities of such a sample, and plot in the top 
panel of Figure [8] the values of T; and derived for each 
HVS. Using equations JT2J and (fH)) , we obtain T = 0.50 and 
<jf. = 0.02. For a a = 5 km s" 1 , we obtain T = 0.50 and of — 
0.03. These errors are comparable with the systematic error 
caused by the flattened disk (see the slight offset of dots from 
the solid line in Fig.[6ji). We have tried different values of the 
gravitational potential parameters (p, q) in the calculations, 
and found that the estimated value of T is always consistent 
with the assumed halo triaxiality (see middle and bottom 
panels in Fig. [5} . 

If the x-axis of the halo potential form an angle r/o (mea- 
sured counter-clockwise) with the reference axis of the ob- 
servations, then (as seen from Fig. [9} the minimum variance 
of Jz,i/J y ,i in a set of simulated samples each in a frame 
at angle rj from the observational reference axis (see eq. I17|l 
occurs for rj = rjo- Our calculations show that an error of 
5° in the estimate of r/o may cause an error of 2% in the 
estimate of T. 

Note that the axis ratios (p, q) of the halo potential are 
degenerate in the defined triaxiality parameter (eq.|5J. After 
determining T, the values of (p, q) could be also obtained by 
using any value of {J x ,J y ,Jz) (e.g., see Fig. |6p), but the 
modeling would be sensitive to halo mass and the shape of 
h alo potential used (s ee also the determination of axis ratios 
m iGnedin et all 12005 by tracing back HVS orbits). 
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Figure 7. The variable (v/r)J z versus (vjr)J y . The solid and dashed lines correspond to J z = —TJ y and J z = (see eq.0, respectively. 
The stars are those shown in Fig. [6] Panel (a) includes only stars far from the plane of the disk, |6* — 90° | > 30°, while Panel (b) includes 
all stars. The effect of the disk potential can be gauged from the spread of the stellar dots from the solid line to the dashed line in Panel 
(b). The magenta and red dots scattered above the dashed line or below the solid line represent bound stars with bent orbits. 
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Figure 8. Simulated triaxiality parameters and their errors from 
a mock sample of HVSs with v > 300 kms -1 , |6» - 90°| > 30°, 
and 55 kpc < r < 200 kpc. An observational error of 3kms _1 is 
assumed for in "measured" one-dimensional velocity. From top to 
bottom, the values of (p,q) used in the calculations are (0.8, 0.7), 
(0.8, 0.6), and (0.8, 0.75), corresponding to a triaxiality T=0.54, 
0.32, and 0.72, respectively (horizontal lines). A statistical anal- 
ysis of the sample using eqs. (1 1 2 I t and (1 1 3 I t yields (f , o-j, )=(0.50, 
0.02), (0.31, 0.01), and (0.64, 0.02), respectively. For clarity, only 
1/3 of the sample points are drawn in the figure. 
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Figure 9. Variance of the triaxiality parameter T" of the simu- 
lated sample (top panel in Fig. [8} in different reference frames as a 
function of rj— rjo- Here rjo is the angle measured counter-clockwise 
from a reference direction (e.g. the line from the Galactic center 
to the Sun) to the a;-axis of the halo potential, and r\ is the angle 
formed by the observational frame with the reference direction. 
The variance has a minimum at r\ — r]o = 0° . 



Finally, it is possible that a few HVSs in the halo 
may be produced by the interactions of stars with an 
IMBH in satellite g alaxie s li ke the Large Magellanic Cloud 
llEdelmann et all 120051 ; iGualandris fc Portegies Zward 
2007), and that these would contaminate the sample ejected 
from the Galactic center. The ejection rates from the 
satellite dwarfs of the MW are expected to be much smaller 
than the rate from Sgr A*, however. Such "satellite" HVSs 
will also have much larger angular momenta in the Galac- 
tocentric rest-frame, and should be easily distinguishable 
from the Galactic center sample. 
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APPENDIX A: CORRECTION OF THE 
EFFECT OF DISK POTENTIAL ON THE 
INFERRED HALO TRIAXIALITY 

As shown in Figure |6j a), the triaxiality obtained through 
T = —Jzjjy has a small offset from the true halo triaxial- 
ity. Part of this offset comes from the effects of the flattened 
disk potential, and part from the fixing of (0, <f>) in the inte- 
gration of equation (|10[) . The offset due to the constant (8, <f>) 
assumption is smaller for stars of higher velocities, as seen 
in Figure |6je) that shows the results after the disk potential 
was removed. Below we provide a method to correct the off- 
set due to the Galactic disk. As in the determination of halo 
triaxiality proposed in this work, even for this correction it 
is not necessary to trace back the orbits of HVSs. 

The contribution to the change in specific angular mo- 
mentum due to the disk potential can be expressed as: 



J, 



y.disk — 



where 

dJy,dis 

df 

dt _ 
dr i 



° djj^disk dt 

dt dr 



dr, 



2 ( d ® disk _ d^disk 
V dz 2 dx 2 



^v 2 + 2$(r ) - 2$(r) : 



(Al) 



(A2) 



(A3) 



r e j is the initial distance of the HVS from the Galactic center 
at ejection, fo and vq are the current position and velocity 
of the HVS, ro = |fo[, and v r is its radial velocity at r. We 
can remove the effect of the disk potential from the variable 
J y by computing J y — J^disk, where 



<Jy , dis 



y , disk 



sm cos cos ( 



GM dis 



(r 2 + «Lk + &disk + 2a disk ^r 2 cos 2 e + b 2 iBk r/ 
dr 



adiskf 



disk 



(A4) 



and where equation ([20} describing the disk potential has 
been used. The angles (9, (f>) are fixed in the integration to 
be present-day values. In addition to the disk potential, we 
also need to assume a halo potential to determine v r , for 
which we may use the spherical part of the halo potential 
by setting (p, q) — (1,1). We have tested this correction and 
found that half of the offset to the true halo triaxiality can 
be corrected. 
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